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ABSTRACT 

Aims. We investigate the environment of the infrared dust bubble S5 1 and search for evidence of 
triggered star formation in its surroundings. 

Methods. We performed a multiwavelength study of the region around S5 1 with data taken from 
large-scale surveys: 2MASS, GLIMPSE, MIPSGAL, IRAS, and MALT90. We analyzed the spec- 
tral profile and the distribution of the molecular gas ( 13 CO, C ls O, HCN, HNC, HCO + , C 2 H, 
N 2 H + , and HC 3 N), and dust in the environment of S51. We used a mid-infrared emission three- 
color image to explore the physical environment and GLIMPSE color-color diagram [5.8]-[8.0] 
versus [3.6]-[4.5] to search for young stellar objects and identify the ionizing star candidates. 
Results. From a comparison of the morphology of the molecular gas and the Spitzer 8.0 //m 
emission, we conclude that the dust bubble is interacting with CO at a kinematic distance of 3.4 
kpc. The bubble S51 structure, carried with shell and front side, is exhibited with 13 CO and C ls O 
emission. Both outflow and inflow may exist in sources in the shell of bubble S51. We discover 
a small bubble G332.646-0.606 (fl,„ = 26", r in = 15", R 0M = 35" and r out = 25") located at the 
northwest border of S51. A water maser, a methanol maser, and IRAS 16158-5055 are located at 
the junction of the two bubbles. Several young stellar objects are distributed along an arc-shaped 
structure near the S51 shell. They may represent a second generation of stars whose formation 
was triggered by the bubble expanding into the molecular gas. 

Key words, infrared: stars — stars: formation — ISM: bubbles — HII regions 

1. Introduction 



Churchwell et al 



(2006 



2007) detected and cataloged about 600 mid-infrared dust (MIR) bubbles 
between longitudes -60° and +60°. The bubbles have bright 8.0 jum shells that enclose bright 24 
yum interiors. The infrared (IR) dust bubbles may be produced by exciting O- and/or B-type stars, 
which are located inside the bubble. The ultraviolet (UV) r adiation from exc iting stars may heat 



2008), which is known 



dust and ionize the gas to form an expanding bubble shell ([Watson et al. 
as the "collect-and-collapse" process. This process can trigger the massive star formation near the 
shell clumps. These bubbles present an important opportunity to study the interaction between HII 
regions and molecular clouds. 
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We selected the IR dust bubble S51 from the catalog of 



Churchwell et al. 



(2006). S51 is a 



complete (clos ed ring) IR dust bubble centered on Z-332.673, fr=-0.618 . It lies in the southern part 



of RCW106 (IMookeriea et al 



2004 



Wong et al. 



2008 



Lo et al 



2009). The radius of the shell is 



1.48' (1.47 pc), and the average thickness of the shell is 0.32' (0.32 pc), while the eccentricity of 
the ellipse is 0.82. The assumed distance is 3.4 kpc. 

MALT90 is a pilot survey conducted with t he Mopra telescope in preparation for the Millimeter 



Astronomy Legacy Team Survey at 90 GHz ([Foster et al 



2011). One of the aims of this multi- 



molecular line mapping in this work is to examine how different molecules correlate with each 
other and with bubble S51. The molecules selected include 13 CO, C 18 Q, HCN, HNC, HCO + , C 2 H, 



N 2 H + , and HC 3 N. li CO has an optical depth of 1 to 12 dWong et alJ 120081) . but C ls O can pro- 



vide optical depth and line profile information (n en , = 3 x 10 3 cm -3 ). MALT90 emission traces 
dense gas, and each molecular transition provides slightly different information. For example, HCN 
(«„,-, = 2 x 10 5 cm~ 3 ) traces high column density and is optically thick; HNC («„•„ = 3 x 10 5 cm -3 ) 
is prevalent in cold gas (n c r g — 2 x 10 5 cw~ 3 ); HCQ+ (n cr j t — A x 1 5 c/n~ 3 ) often shows infall sig- 
natures and outflow wings ( Rawlings et alj|2004 ; Fuller etal]|2005 ): N 2 H + (n crit = 4x 10 5 cm -3 ) is 



more resistant to freeze-out onto grains than the carbon-bearing species (Berginetal EoOll) : HC 3 N 
is a good tracer of hot core chemistry. 

In this work, we report a multiwavelength study of the environment surrounding the IR dust 
bubble S5 1 . We aim to explore its surrounding interstellar medium (ISM) and search for signatures 
of star formation. We describe the data used in Sectj2] the results and discussion are presented in 
SectfSJ Sect|4]summarizes the results. 



2. Data 



2003 



MALT90 (Foster et al. 



2011) 



Churchwell et al. 2009), MIPSGAL (Carey et al. 2009), IRAS (Neugebaueret al 



We analyzed IR and millimeter wavelength da ta extracted from the f ollowing larg e-scale surveys: 
the Two Micron All Sky Su rvey (2MASs|] (S krutskie et al. 2006 ) GLIMPSE (Beniamin et al 



1984), and 



GLIMPSE is a MIR survey of the inner Galaxy performed with the Spitzer Space Telescope. We 
used the mosaicked images from GLIMPSE and the GLIMPSE Point-Source Catalog (GPSC) in 
the Spitzer-IRAC ( 3.6, 4.5, 5.8 and 8.0 yum). IRAC has an angular resolution between 1.5" and 1.9" 



(Fazi o et al. 



2004 



Werner et al. 



2004). MIPSGAL is a survey of the same region as GLIMPSE, 
using the MIPS instrument (24 and 70 jum) on Spitzer. The MIPSGAL resolution is 6" at 24 jim. 

The Mopra 22-m radio telescope has a full width at half maximum (FWHM) beam size of ~ 33" 
at ~ 1 10 GHz for 13 CO 7=1-0 and C 18 0/= 1-0 transitions. The correlator was configured with 



per channel over a useable velocity bandwidth of ~120 km s 1 (1L0 et al. 



Lo et al. 


2009; 


Bains et al. 


2006 



Wong et al.ll2008l) . The observing was set up so that the central channel corresponded to ~50 km 
s _1 , at which the velocity of the emission from the GMC complex is centered. 



1 2MASS is a joint project of the University of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and 
the National Science Foundation. 
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332.7 332.65 332.7 332.65 332.7 332.65 332.7 332.65 332.7 332.65 

Galactic longitude 

Fig. 1. Integrated velocity contours of the 13 CO emission every 1.0 km s superimposed on the 
GLIMPSE 8.0 /urn grayscale. The contours are at multiples of the 13% level of each of the emission 
peaks. The peak of each mosaic from -62 to -43 km s _1 is 0.718, 1.632, 3.584, 7.503, 9.191, 
11.388, 14.901, 16.428, 10.915, 9.673, 12.091, 15.416, 16.916, 16.467, 14.378, 12.772, 6.043, 



3.018, 1.535, and 2.631 K km s , respectively. The green symbols "a'' 



' and indicate the 



positions of the water maser, the methanol maser, and IRAS 16158-5055, respectively. 



The on-the-fly (OTF) mapping mode of Mopra was used for MALT90. Maps were made with 
the beam center running on a 3.4' x 3.4' grid. The scan rate was 3.92" per second. The map was 
made with 12" spacing between adjacent rows, giving 17 rows per map. Because the Mopra beam 
at 90 GHz is 38", this row spacing provides redundancy in the map. The full 8 GHz bandwidth 
of the Mopra spectrometer (MOPS) was split into 16 zoom bands of 138 MHz each, providing 
a velocity resolution of ~0.11 km s in each band. In this work, six MALT90 pilot survey lines 
(HCN(7 = 1-0), HNC (7 = 1-0),HCO + (7 = 1-0), C 2 H (/V = 1 -0, 7 = 3/2- 1/2,F = 2 - 1), 
N 2 H + (7=1-0) and HC 3 N (7 = 10-9)) were used to trace the environment of bubble S51. 

For the MOPRA observations, OFF positions were chosen at ±1° in Galactic latitude away from 
the plane (positive offset for sources at positive Galactic latitude and vice versa), and a single OFF 
position was observed for every two scan rows. Pointing on SiO masers was performed every 1- 
1.5 hour, maintaining pointing precision to better than about 10". The brightness temperature Tmb 
is related to the antenna temp erature 7! by Tm b — J)Jt\v, where n v is the frequency-dependent 



Ladd etal 



(20051) and 



Lo et al 



(2009), the main beam efficiency at 



beam efficiency. According to '. 
86 GHz is ?786GHz = 0.49, and at 110 GHz is jjwochz - 0.44. The results presented in this paper 
are i n terms of main bright temperature Tmb(K). The I3 CO and C I8 data cubes were obtained 
( 2009 ). MALT90 data cubes were downloaded from the online archivj^. The 13 CO, 



from 



Lo et al. 



http://atoa.atnf.csiro.au/MALT90/ 
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332.7 332.65 332.7 332.65 332.7 332.65 332.7 332.65 332.7 332.65 

Galactic longitude 

Fig. 2. Integrated velocity contours of the C 18 emission every 1.0 km s superimposed on the 
GLIMPSE 8.0 yi/m grayscale. The contours are at multiples of the 13% level of each of the emission 
peaks. The peak of each mosaic from -62 to -43 km s" 1 is 0.518, 0.590, 0.870, 2.625, 4.342, 4.081, 
3.663, 2.814, 1.627, 1.529, 2.671, 3.697, 4.914, 4.572, 4.514, 3.756, 1.524, 0.575, 0.558, and 0.685 



K km s , respectively. The green symbols "a'' 



and indicate the positions of the water 



maser, the methanol maser, and IRAS 16158-5055, respectively. 



C lg O and MALT90 data cubes were processed with CLASS and GREG in the GILDAS software 
packagj^. 



3. Results and discussion 

The clue of this section is as follows. Firstly, by analyzing the channel maps (FigsQ] and |2| of 
bubble S51 with 13 CO and C lg O emission, we found that bubble S51 has shell and front side 
characters strongly associated with 8.0 fim emission. The integration intensity maps of [-62.0 - 
43.], [-62.0 -53.0] and [-53.0 -43.0] components are shown in Fig[3] In addition, to examine the 
clumps of bubble S51, we used MALT90 multimolecular line mapping (Fig lA.ll in Appendix lAl 
to compare them with each other and with IR emission distribution. The molecules include HCN, 
HNC, HCO + , C2H, N 2 H + , and HC3N. And we also show the spectra of these molecules in Fig]4] 
These spectra indicate that there are signs of inflow or outflow, whose velocity channels are shown 
in Fig|5] Furthermore, we show two three-color images of MIR emission of IR dust bubble S51 
with 3.6-4.5-8.0 fim and 4.5-8.0-24 fim in Fig|6] so that we can understand the IR structure of 
bubble S51. Finally, we used GLIMPSE color-color diagram [5.8]-[8.0] versus [3.6]-[4.5] (FigHl 
and spectral energy distribution (SED) fitting (FigJS) to search for young stars and identify exciting 
star candidates, whose distribution is shown in Fig|9] 
3 http://iram.fr/IRAMFR/GILDAS/ 
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332.7 332.65 332.7 332.65 332.7 332.65 
Galactic longitude 



Fig. 3. Integrated velocity contours of the 13 CO (Upper three mosaics) and C ls O (Lower three 
mosaics) emission superimposed on the GLIMPSE 8.0 //m grayscale. The contours are at multiples 
of the 9% and 12% level of each of the 13 CO and C ls O emission peaks, respectively. The peaks of 
upper three mosaics from left to right are 123.136, 59.055, and 93.380 K km s _1 ; the peaks of the 
lower three mosaics from left to right are 25.598, 18.456, and 20.458 K km s _1 , respectively. The 
integrated velocity range is indicated in each panel. The green symbols "a", "■" and "it" indicate 
the positions of the water maser, the methanol maser, and IRAS 16158-5055, respectively. 

3. 1 . The shell and front side of bubble S51 

We inspected the molecular gas around S5 1 from the Mopra 13 CO and C I8 emission in the whole 
velocity range and found two independent velocity components ([-62.0 -53.0] and [-53.0 -43.0] 
km s" 1 ) along the line of sight. FigsQ] and [2] show the integrated velocity maps of the 13 CO and 
C ls O emission every 1.0 km s _1 between -62.5 and -42.5 km s _1 . The 8.0 fim grayscale shows the 
position and size of S51. 

Comparing the molecular and MIR morphologies, there are many interesting places among the 
[-62.0 -53.0], [-53.0 -43.0] km s _1 components and 8.0 fim emission, shown in Figs[T]|2]and[3] The 
morphology and peaks of [-53.0 -43.0] km s _1 component clearly correlate with 8.0 fim emission. 
Accordingly, the [-53.0 -43.0] km s _I component may be the shell of bubble S51. On the other 
hand, looking at the -58.0 km s _1 panel in Figj2] we find that the contours fit quite well within 
the 8.0 fim emission ring. The integration velocity map of the [-62.0 -53.0] km s _1 component in 
FigEJalso shows that there is a good correlation between C 18 contours and 8.0 fim emission. We 
suggest that this cloud could be blueshifted gas associated with the front side of bubble S5 1 . Since 
the -48.0 km s _1 cloud is well associated with bubble shell, this bubble would be expanding at 
roughly 10.0 km s _1 along the line of sight. However, there is no redshifted counterpart. 

From the 13 CO and C 18 spectra of Figj4] we can also see that there are two velocity compo- 
nents at the locations of sources A, B and C (see Fig|9]l, where source A is the core of eastern lobe 
of S51, source B lies at the front side or inside S51, and source C is the core of western lobe of 
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Fig. 4. MALT90 spectra of sources A, B and C, whose positions are represented in Fig|9]wifh "A" 
(7=332.693, fc=-0.610),"B" (7=332.677, £=-0.619) and "C" (Z=332.645, £=-0.608 ), respectively. 
The red dashed lines indicate the peak velocities of the C I8 emission. The green curves are Gauss- 
fit lines. HCN and N2H + have hyperfine splitting, shown as a black bar. 

S51. It is obvious that the spectral intensity of the [-62.0 -53.0] km s _1 component at the front side 
of bubble is stronger than that of [-53.0 -43.0] km s _1 . At the shell of bubble S51, however, it is 
just the reverse. This case is also evidence that bubble S51 has a shell and front side structure. 

For the [-53.0 -43.0] km s _1 component, the optical thin C is O emission contours and spectra 
(Figf2]and|4]i s how that th e systematic velocity is — 50.0 km s _1 . Assuming the galactic rotation 



model of Fichetal 



(1989) ( with j?p = 8.5 kpc and = 220 km s" 1 ), we obtain kinematic distances 



of either 3 .4 or 1 1 .7 kpc. lLockmanl (1 19791) suggested that RC W 1 06 lies the near kinematic distance 
after examining H2CO absorption spectra seen against the HII region continuum. Therefore, we 
adopted the near kinematic distance of 3.4 kpc in this work. 



3.2. Opacity, excitation temperature, and column density 



We used the approach of IWong et al.1 (120081) to derive the opacity, excitation temperature, and 
column density of 13 CO and C 18 within the box shown in Fig. lA.ll of Appendix lAl The size of 
the box is 4.3' x 4.3' (where 4.3' is the diameter of S51). The velocity range included is just V - 
—53.0 to —43.0 km s" 1 of bubble shell, and the box is centered on Z=332.673, £=-0.618. 

The relation be tween opacities and the ratio of 13 CO to C 18 main-beam brightness temperature 



(Myers et al. 1983) is 



Tmb( 13 CO) 



1 - exp(-Ti 3 ) 1 - exp{-lATii) 



T MB (C l *0) 1 - exp(-Tis) 1 - exp(-Tis) 



(1) 
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Equation (Q} assumes T13 = 7.4ti8 (IWilson & Rood 



1994 



Won g et al. 



2008). Furthermore, the 



excitation temperature T ex is derived from the equation of radiative transfer: 

| T MB = f[J(T ex ) - J(T hg )][\ - exp(-T)] 
\ J(T) = T /[exp(T /T) - I] 



(2) 



where 



13 CO (IWong et al.ll2008l) . And then, assumin 



ing factor, T/, g = 2.7 K and Tq = hv/k = 5.29 K for the 1-0 transition of 



tical thin, we obtain the molecular 13 CO 



column density A^( 13 CO) from the relation (Bou rke et alll997l) 



NCCO),, 



Tex + 0.88 



2.42 x 10 1 



T MB ( u co)dv. 



(3) 



1 - exp(.-5.29/T ex ) J(T ex ) - J(T bg ) 

Finally, assuming [H2/ 13 CO] abundance ratio is 7xl0 5 dFrerking et al.l 19821) . the molecular hydro 
gen column NiHi) was calculated. The molecular cloud mass was estimated at M S heU ~ 1.3xlO 4 M 
for the bubble shell at a distance of Dssi ~ 3.4 kpc. 



3.3. Clump analysis 

To examine the clumps of bubble S51, we used multimolecular line mappings to compare them 
with each other and with the IR emission distribution. 

In the box shown in Fig. IA.1I the maximum and average column densities N( I3 CO) are 
9.72xl0 17 and 5.98xl0 16 cm -2 , respectively. This figure also shows the main beam brightness tem- 
perature Tmb, the brightness temperature ratio Tmb( 13 CO)/Tmb(C 1& 0), the optical depth t( 13 CO), 
and the excitation temperature r ev ( 13 CO). The distribution of these parameters (shown in color) 
are compared with contours showing the integrated intensity of different molecular transitions. 
The average r MB ( I3 CO)/r MB (C 18 0), t( 13 CO) and r„( 13 CO) are 4.49, 1.76 and 20.91, respec- 
tively. The peak positions of I3 CO, C 18 and 8.0 pm emission are spatially coincident, and are 
situated at the east and west edge of the shell. The molecular clouds are extended around the 8.0 
pm emission peaks. At the south border of bubble S51, we can also see the shell of the filament 
structure from 8.0 pm emission in Fig(6] The 13 CO and C 18 emission contours are distributed 
along the shell of the filament structure in Figf5] We also see from Fig lA.ll that the distribution of 
r MB ( 13 CO)/r M i;(C 18 0), t( 18 CO), 7/ eA ( 13 CO) and N( 13 CO) shows an arc-like morphology around 
the north edge of S5 1 . However, it should be noted that these derived properties are not independent 
of each other. 

HCN, HNC, HCO + , C 2 H, and N 2 H + trace the position of highest density, while HC 3 N traces 



the location of the hot core (IFoster et al 



201 lb . We chose three sources (A, B and C) to analyze 



in Fig|9l their spectra are shown in Figj4] Fig JA. II shows that a potentially dense and hot core is 
mainly located in the eastern and western shells. The spectra shown in Fig@]for the [-53.0 -43.0] 
km s" 1 component demonstrate that source B inside S51 has a lower brightness temperature than 
sources A and C in the S51 shell. At source B, however, the MALT90 spectra of the [-53.0 -43.0] 
km s -1 component show the blueshifted velocity [-62.0 -53.0] km s~\ which is evidence that this 
is the front side of bubble S51. This suggests that the extended cloud of the [-53.0 -43.0] km s~' 
component is scarce inside bubble S51, and the front side of bubble S51 has a rich abundance. 
However, the back side of bubble S5 1 is not shown in the MALT90 spectra. They also show that 
bubble S51 might be expanding at 10 km s~' compared to the bubble shell along the line of sight. 
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332.68 332.66 332.64 332.62 
Galactic longitude 




332.68 332.66 332.64 332.62 
Galactic longitude 



Fig. 5. Integrated velocity contours of HNC (top panel) and HCO + (bottom panel) emission super- 
imposed on the GLIMPSE 8.0 /im grayscale. The contours are at multiples of 11%, 11%, 12%, 
and 12% level of the blue lobe of HNC, the red lobe of HNC, the blue lobe of HCO + , and the red 
lobe ofHCO + emission peaks, which are 9.539, 8.618, 12.861, and 10.800 K km s" 1 , respectively. 
Source C is indicated with "C" in the image. The symbols "a", "■" and "jt" indicate the positions 
of the water maser, the methanol maser, and IRAS 16158-5055, respectively. 



In addition, we note that IRAS 16158-5055 is closer to the peaks of N 2 H + and HC3N than the water 
maser and the methanol maser. 
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Fig. 6. MIR emission of IR dust bubble S5 1 . Top: Spitzer-IRAC three-color image (3.6 fim = red, 
4.5 fim = blue and 8.0 pm = green); Bottom: Spitzer-IRAC and Spitzer-MIPSGAL three-color 
image (4.5 /mi = blue, 8.0 pm = green and 24 pm = red). Note that the 24 pm emission is saturated 
at the center of image. 



3.4. Inflow or outflow character 

As shown in FigO there is much trace of active star formation located in the shell of S51. The 
western shell of S51 contains a UCHII region, a water maser (/=332.653, £=- 0.621), a methanol 



maser (Z=332.660 , £=-0.618), and IRAS 16158-5055 (Z=332.653, £=-0.612) (Wal sh et al. 



1997 



Breenetal 



20071) . The water maser spectrum has three peaks, with fluxes of ~1 Jy for the -58.5 
km s -1 component, ~6 Jy for -47.5 km s~', and ~30 Jy for -45.2 km s~', respectively (Wals h et al 
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Fig. 7. GLIMPSE color-color diagram [5.8]-[8.0] versus [3.6]-[4.5] for sources in and around S51. 
The symbols "a", "•" and "■" indicate the different stars of class I, II, III and others, respec- 
tively. The classifi cation of class I, II and III indicates the stellar evolutionary stage as defined by 



Allen et al 



(2004) 



1997ft . T he methanol maser spectrum has a velocity range of -52 ~ -50 km s _1 , and the peak flux 
is 7.1 Jy (iBreen et al.ll2007l) . This means that the [-53.0 -43.0] km s _1 component of the S51 and 
not the [-62.0 -53.0] km s _1 component of background is strongly correlated with the methanol and 
water masers. 

On the west side of S51, the position of source C is marked in Fig|5] Fig0] shows the 13 CO, 
C ls O and MALT90 spectra of source C. We marked the splittings of the hyperfine structure of 
HCN and N2H + with several black bars. The [-62.0 -53.0] km s" 1 component is very faint in 13 CO, 
C 18 emission, as well as HCN, HNC, HCO + , C 2 H, and HC 3 N. Consequently the [-62.0 -53.0] 
km s" 1 component seems unrelated to the [-53.0 -43.0] km s _1 component in this cluster. Fig0]also 
shows that the optically thick spectra of HCN, HNC and HCO + lines have absorption dips at the 
peak velocity (-49.5 km s~') traced by C ls O, C 2 H, N 2 H + and HC 3 N. The blue profile is stronger 
than the red profile shown from the HC N spectrum. Such a pro file can be produced by inflow with 



cold gas wings and a hot central core (Ma rdones et al 



1997). Moreover, HNC and HCO + have 



broad line wings corresponding to N 2 H + and HC3N, therefore source C could also be produced by 
outflow. Fig|5] shows the emission in two different velocity channels. The velocities of the water 
maser and the methanol maser are similar to the more negative velocity gas. We speculate that 
since this region is located at the junction of bubble S51 and bubble G332. 646-0. 606 (referred to 
Sect l3.5b . it is plausible that the interactions between the two bubbles or between the HII region 
and the molecular cloud have triggered this phenomenon. 



3.5. IR structure of S51 



Fig|6] shows two Spitzer three-color images of S51. Both figures clearly show the PDR visible in 
8.0 fim (in green) emission, which originates mainly in the PAHs. Since these large molecules 
are destroyed inside the ionized region, the PAH emission delineates the HII region boundaries, 
and the molecular c l ouds are excited in the PDR by the radiation leaking from the HII region 



(IPetriella et al 



2010; 



Pomares et al 



2009). The 24 /jm emission (bottom panel in Fig|6]l appears 
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Fig. 8. SED fits of exciting star source 10 (7=332.6664, b=-0.6266). The flux density of fitted 
data are the 2Mass H-band 5.01(0.22), the K-band 17.81(0.44), IRAC 3.6 30.16(1.26), 4.5 jum 
20.93(1.48), 5.8 fim 25.12(1.53), and 8.0 fim 22.15(3.46) mJy, respectively. 



to be inside bubble S51, and corresponds to hot dust. It is likely that O- and/or early B-type stars 
produced the bubble shell of this HII region, with hot dust located inside the bubble. The 3.6 //m 
emission (in red, top panel in Figj6j shows the positions of the brightest stars. 

We also note a small bubble G332. 646-0. 606, which is centered on /=332.646, £>=-0.606 and 
located at the northwest edge of S51. The semimajor (Ri„) and semiminor (r,„) axes of the inner el- 
lipse are 26 and 15 arcsec, respectively; the semimajor (R oul ) and semiminor (r out ) axes of the outer 
ellipse are 35 and 25 arcsec, respectively. From 8.0 //m emission of bubble G332. 646-0. 606, the 
elliptical morphologic PDR shows a sketch of the bubble. Note that the 24 fim emission is saturated 
inside the small bubble G332. 646-0.606. At the junction between bubble S5 1 and G332. 646-0.606, 
the IR emission is brighter than the rest of the shell; this is also true for the I3 CO, C 18 and 
MALT90 emission lines. IRAS 16158-5055 is also located at this position. It seems likely that this 
junction region is the birth space of several YSOs. In addition, the molecular gas of the bubble shell 
exhibits several clumps along the PDR (FigO. The distribution and morphology of this material 
suggests that a collect-and-collapse process may be occurring. 



3.6. Search for young stars and identifying exciting stars 

The distribution of the IR point sources in the surroundings of S5 1 provides some signs of star 
formation. Fig|7] shows the [5.8]-[8.0] versus [3.6]-[4.5] color-color diagram for the sources ex- 
tracted from the GLIMPSE Point Source Catalog in the Spitzer-IRAC bands in and around S5 1 . 
We only considered sources with detection in the four Spitzer-IRAC bands. In Fig(7J class I stars 
are protostars with circumstellar envelopes; class II are disk-dominated objects; and class III are 
m ain sequence and g iant stars. The classification indicates the stellar evolutionary stage as defined 
by 



Allen et al 



(12004 . 

Class I and class II stars are chosen to be the YSO candidates. The YSO candidates are drawn 
with symbols "x" and "a" in Fig|9] among which sources 4, 7, 8, 9, 13 and 17 near the shell of 
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332.7 332.65 
Galactic longitude 



Fig. 9. Integrated velocity contours of the 18 CO emission superimposed on the GLIMPSE 3.6 fim 
grayscale. The contours are at multiples of 12% level of the C 18 emission peaks 93.380 K km s _1 , 
and the integrated velocity range is from -53.0 to -43.0 km s _1 . Two dashed ellipses show the PAH 
morphology of the two bubbles. Letters "A" (7=332.693, fc=-0.610),"B" (7=332.677, £=-0.619) and 
"C" (7=332.645, £=-0.608) indicate the chosen positions of three sources. The class I, II and III 
stars are indicated with symbols "x", "a" and "O". The symbols "a", "■" and indicate the 
positions of the water maser, the methanol maser, and IRAS point source, respectively. 



bubble are made up of an arc-shaped distribution. This may be caused by interaction between the 
HII region and molecular clouds, and be the production of a collect-and-collapse process. 

We also tried to search for and identify exciting stars of S51 . However, only sources 10 and 14 
are located within the bubble. Source 14 belongs to a class I star, therefore we did not consider it to 
be an exciting star. Source 10 is a cl ass II star, therefore we identified it in more detail with an SED 
fitting, using the tool developed by Robitaille et al. (2007), which is available onlineQ. The good 
fitting models were selected according to the condition v 2 - v 2 . < 3, where v 2 . is the minimum 

to to /l A mm ' A-min 

value of^- 2 among all models. 

In Figj8] we show the results f or the SEP fitting of our source 10 to the fluxes obtained from the 
GLIMPSE Point Sour ce Catalog ( Horaet alj|2008 ) and the 2MASS All-Sky Point Source Catalog 
(ISkrutskie et al.l 120061) . We fit these sources allowing the extinction to range from to 30 mag 
and the distance to range from 3.0 to 4.0 kpc. The SED output shows that source 10 in class II 
is Menv/M* < 10~ 6 yr~' and M disk [Mi, < 10~ 6 ; the temperature of th e center star is 31000 K, i n 
agreement with the effective temperature expected for an O-type star ( Schaerer & de Koteilll997 ). 
Therefore, source 10 may be a candidate of the ionization of S51. Stellar winds emit from exciting 



http://caravan.astro.wisc.edu/protostars/ 
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stars. In addition, the presence of a central cavity in the distribution of the ionized gas may be a 
signature of the activities of stellar winds. The location of source 10 inside this cavity provides 
additional support for the assumption that this star is the most likely candidate to have created S5 1 . 
Spectroscopic observation is needed to confirm this. 

4. Summary 

We have investigated the environment of the IR dust bubble S51 with several spectra ( 13 CO, C ls O, 
HCN, HNC, HCO + , C 2 H, N 2 H + , and HC 3 N) and IR emission. The main results can be summarized 
as follows. 

(1) We have distinguished two independent velocity components associated with S51 along the 
line of sight. One component belongs to the shell of bubble S51, and another may be the front side 
of bubble S51. 

(2) We suggest that source C (7=332.645, £>=-0.608) shows evidence of either outflow or inflow, 
located on the northwest border of S5 1 . 

(3) Next to the west border of bubble S51, we found a small bubble G332.646-0.606, whose 
dimensions are = 26", r in = 15", R out = 35" and r out = 25". 

(4) MALT90 emission contours appear to be correlated with the S5 1 shell. The water maser, 
the methanol maser, and IRAS 16158-5055 are located in the shell. The YSO distribution is also 
correlated with 13 CO, C ls O and 8.0 yum emission. 
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Appendix A: Molecular line maps of bubble S51 
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Fig. A.l. Integrated velocity contours of MALT90 emissions (HCN, HNC, HCO + , C 2 H, N 2 H + , 
and HC3N) superimposed on six color images of different parameter distributions (Tmb( 13 CO), 
ratio = r MB ( 13 CO)/r MB (C 18 0), r( 18 CO), r eA ( 13 CO), N( 13 CO)[xl0 16 ] and 8 /an). The contours are 
at multiples of the 12%, 13%, 12%, 16%, 14%, and 15% level of each of the HCN, HNC, HCO + , 
C 2 H, N 2 H + , and HC 3 N emission peaks, which are 9.573, 7.953, 9.781, 4.791, 10.490, and 2.555 K 
km s , respectively. The integrated velocity range is from — 53.0 to — 43.0 km s for all contours. 
The black symbols "a", "■" and indicate the positions of the water maser, the methanol maser, 
and IRAS 16158-5055, respectively. 



